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The main feature in the elastic neutron scattering of La2- a; Sr :E Cu04 is the existence of incommen- 
surate peaks with positions that jump from 45° to 0° at 5% doping. We show that the spiral state 
of the t — t' — t'—J model with realistic parameters describes this data perfectly. We explain why 
in the insulator the peak is at 45° while it switches to 0° precisely at the insulator-metal transition. 
The calculated positions of the peaks are in agreement with the data in both phases. 
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Introduction. The phase diagram of L^-zSr^CuC^ 
(LSCO) shows that the magnetic state changes dramat- 
ically with Sr doping. The parent compound La 2 Cu04 
exhibits three-dimensional long-range antifcrromagnetic 
(AF) order below 325K The Neel order disappears 
at Sr concentration x ps 0.02, however two-dimensional 
(2D) short-range AF correlations exist at any doping 
(see also Ref. yj for a review). At x < 0.055 the system 
exhibits only hopping conductivity and behaves like an 
Anderson insulator, while the usual dc conductivity as 
well as superconductivity appear at x > 0.055 0,0. 

Static magnetic ordering at very low temperatures has 
been observed both for x < 0.055 and x > 0.055. The 
elastic neutron scattering peak is close to the AF position 
Qo = (7r, 71"), but is shifted from this position by SQ: 
Q = Qo + SQ. We set the lattice spacing a = 1. This 
shift indicates a one-dimensional incommensurate spin 
modulation. The dependence of the shift on doping has 
been studied in the superconducting phase Q, as well 
as in the insulating phase H, H, • These studies have 
revealed the following remarkably simple dependence of 
the elastic peak shift on doping x (see Fig. 7 in 0): 

0.055 < x < 0.12 : SQ ps 2x{±tt, 0) or SQ ps 2x(0, ±tt) , 
0.02 < x < 0.055 : SQ ps V2x(±tt, =pr) . (1) 

Thus the ID incommensurate spin modulation is propor- 
tional to doping and the direction jumps from 45° to 0° 
exactly at the point of the insulator-metal transition. 

One of the early proposals made by Shraiman and Sig- 
gia in Ref. |8j, and later explored in the context of the 
Hubbard and the t-J models @, ES El El El 13 EH , 
was that for small doping the collinear Neel order gives 
way to a non-collinear spiral state. There is a gain in 
energy since the holes can hop easier in a spiral back- 
ground. However the issue of stability of the spiral state 
remained rather controversial. Using chiral perturbation 
theory E3 we have recently revisited the problem of sta- 
bility of the spiral state in the extended t — t' — t" — J 
model jlFij . and have found that the uniform (1,0) spiral 



state is stable (at low doping) above some critical val- 
ues of t' , t" . The stability is due to quantum fluctuations 
(order from disorder effect). Even more importantly, su- 
perconductivity coexists with the spiral order. The start- 
ing point of the approach EH is the ground state of the 
Heisenberg model which incorporates all spin quantum 
fluctuations. The chiral perturbation theory allows a reg- 
ular calculation of all physical quantities in the leading 
order approximation in powers of doping x. Sublead- 
ing powers of x depend on the short-range dynamics and 
hence cannot be calculated without uncontrolled approx- 
imations. Therefore the approach is parametrically well 
justified in the limit x <C 1. The phase diagram of the 
t — t' — t" — J model obtained in Refs. E3 is presented 
in Fig. 1. From the Raman data 17] J ss 125meV and 
we set t/ J = 3.1, following the calculations of Andersen 
et al E3- The values t' ps —0.5 J, t" ps 0.3 J for LSCO 
and t' ps -0.8J, t" ps 0.6J for YBCO are taken from the 
same calculation. From now on we measure all energies 
in units of J (J = 1). The matrix elements t' and t" 
are small compared to t, but nevertheless are crucially 
important for the stability because they influence sub- 
stantially the hole dispersion. The pitch of the uniform 
superconducting spiral state is 

Zf 

5Q= — x (1,0) « 5.8 x (1,0) , (2) 

Ps 

where p s ps 0.18 is the spin stiffness of the Heisenberg 
model and Z ps 0.34 is the quasiparticle residue. The 
residue de pen ds weakly on t' and t", and 0.34 is the value 
for LSCO E3- Notice that Eq. © is in very good agree- 
ment with the data in the metallic phase, x > 0.055. 
However different analysis is needed to explain the data 
JTJ in the insulating region, x < 0.055. 

In what follows we will consider mostly the insulating 
phase. The possibility of spiral ordering in the insula- 
tor has been stressed recently by Hasselmann et al E3- 
However, the dynamical origin of the spiral, the pitch 
of the spiral as well as the "jump" of the spiral direc- 
tion at the insulator-metal transition still remain unex- 
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FIG. 1: The phase diagram of the t—t'—t" — J model at t/J — 
3.1 and small uniform doping 15]. Points corresponding to 
LSCO and YBCO are shown. The line separates the stable 
spiral phase and the unstable region (labeled "?"). 

plained. It is the purpose of the present work to explain 
all these phenomena. In essence our idea is the follow- 
ing. The spiral (l,0)-state of the conductor has lower 
energy than that of the (l,l)-state only due to the pres- 
ence of the Fermi motion energy . On the other hand 
without the Fermi motion contribution the (l,l)-energy 
is lower than the (l,0)-energy. In the region x < 0.055 
each hole is localized near its Sr ion, hence there is no 
Pauli blocking (Fermi energy) and the system immedi- 
ately crosses over to the (1,1) spiral state. Following 
this line of reasoning we will demonstrate below that the 
whole variety of experimental data in the insulator and 
across the insulator-metal boundary can be consistently 
explained. In the insulator the ground state is strongly 
non-uniform, with the holes trapped in the vicinity of 
the randomly distributed Sr ions, as confirmed by the 
variable range hopping (V RH) behavior of the dc con- 
ductivity for x < 0.055 |3J, |20j. Due to the presence of 
the Coulomb potential (see below) and disorder, it seems 
likely that the transition to the metallic (uniform) phase 
at x — 0.055 is of the density-driven percolation type. 
Thus we adopt this point of view although the exact na- 
ture of the transition is not crucially important for our 
analysis. 

Coulomb trapping of holes. Let us consider first a 
single Sr ion with a single hole in an AF background. 
The hole is trapped near Sr by the Coulomb poten- 
tial e 2 /(e e vV 2 + d 2 ) ki e 2 /(e e r) where d is the distance 
from the CuC>2 plane to the Sr ion and e e is the ef- 
fective dielectric constant. For zero doping e e ~ 30 
and it increases with doping, as discussed in [3j. In 
momentum space the hole is localized near one of the 
points ko = (±7r/2, ±7r/2), which are the centers of 
the four faces of the magnetic Brillouin zone (MBZ). In 
the vicinity of these points the dispersion is quadratic: 
£k ~ ^-kf + ^ffc 2 , wnere k is defined with respect to ko, 
and ki is perpendicular to the face of the MBZ, while &2 
is parallel to it. For values of t' and t" corresponding to 



LSCO we find that the dispersion is practically isotropic 
Pi ~ Pi = P 2.2 [l^. Since the lattice spacing is about 
3.85A this value corresponds to an effective mass of about 
two free electron masses, in agreement with the optical 
conductivity data The solution of the Schrodinger 
equation 

(-/3V 2 /2 - e 2 /(e e r)) x = e X (3) 

determines the ground state wave function and the 
ground state energy of the localized hole: 

X(r) = y/2/7TKe- Kr , e = -/3k 2 /2 (4) 

where K — 2e 2 /(e e /3). Using the hopping conductivity 
data at very low doping (x — 0.002) [21| we estimate 
the inverse size k « 0.4. At higher doping (but still 
in the insulator) the value of k might decrease slightly. 
Thus throughout the insulating phase k is small and this 
justifies the semiclassical approximation we use below. 
Note that the semiclassical approximation is used only 
with respect to the kinetic energy of the hole but not for 
the spin (we do not use the 1 /S expansion and account for 
all spin quantum fluctuations via the chiral perturbation 
theory [poT]'). 

Spiral induced by a single trapped hole. Equations @ 
and Q assume a rigid antiferromagnetic background. 
However one can gain energy relaxing the background 
into the spiral state. In the spiral state there are still 
two sublattices, sublattice "up" and sublattice "down", 
but the spin at every site of each sublattice is rotated by 
an angle Qi with respect to the orientation at r = oo 

|i) = e «('t)"»-'/2| t), i S "up" sublattice, 

|j) = e ^(r J ) m - t r/2| ^ j g " down " sub i att i C e. (5) 

Here m = (cos a, sin a, 0) with arbitrary a is the "di- 
rector" of the spiral which is orthogonal to the magne- 
tization plane. Note that directions in spin space are 
completely independent of directions in coordinate space. 
The wave function of the hole ip(r) has two components 
corresponding to up and down sublattices. The total en- 
ergy is of the form 

E = jd*r{^{V9f+ (6) 

, / _#vi_iL y/2Zte- ta (e ■ V0)\ , 1 

*'<"(A.^ ')«"} • 

where e = (-i=,±-i=) is a unit vector orthogonal to a 
given face of the MBZ. Let us search for a solution in the 
form 

V»(r) = -L (_* ia ) X (r) (7) 

where x( r ) is given by Eq. Variation of the energy 
© with respect to 9 leads to the following equation 

Zf 

V 2 = V2— (e-V) X 2 (r) . (8) 

Ps 
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The solution of |[SJ is 

Zt (e-r) 



reads 



\f2-Kps 



[l- e- 2Kr (l + 2nr)] 



(9) 



Substitution of this solution together with Q) and in 
Eq. JBJ) gives the following total energy 



E 



(3 Z 2 t< , , 

77 - 1 « - 2 — « 

2 47rp s ) e„ 



(10) 



Minimizing this energy one finds k = — — 
Z 2 t 2 /(2np s )]. However we do not use directly this ex- 
pression since the effective dielectric constant e e is not 
known accurately enough. Instead we rely on estimates 
for k which directly follow from the hopping conductiv- 
ity as discussed after Eq. As can be easily seen from 
Eqs. (|5I9(I . at distances r <C 1/k our solution describes 
a (1,±1) spiral, while in the opposite limit an effective 
dipole is formed (see discussion below). The solution is 
a variational one because we have used the ansatz Q). 
Even though one can easily derive an exact equation for 
X which can be solved numerically, this is not necessary 
since for our purposes the details of the charge distribu- 
tion are not important. 

We emphasize that the Coulomb trapping of the hole 
is crucially important. Without such trapping the hole 
is delocalized and a single delocalized hole does not gen- 
erate a static spiral. This is qualitatively different from 
the arguments of Ref. 22] . The solution 10 does not 
carry any topological numbers, and consequently, un- 
like the model used in Ref. |23j, our solution is not a 
skyrmion. Other topological reasons for "self-trapping" 
of holes have also been given |24j. however we pursue 
the Coulomb trapping picture since it unambiguously fol- 
lows from the parametrically justified analysis of the t-J 
model. 

The solution 10),® depends on the spiral "director" 
which is a purely classical variable and the energy is in- 
dependent of it. It is unlikely that a finite system has an 
exactly degenerate ground state (spontaneous violation 
of symmetry). This means that higher orders in k (k 4 - 
corrections to the semiclassical solution) may give rise to 
a kinetic energy for the director m and hence to quantum 
rotations of m, lifting the degeneracy. Another quantum 
effect is tunneling from one pocket in momentum space to 
another. However, the quantum corrections are not im- 
portant for understanding the properties of LSCO since 
at finite concentration of impurities the interaction be- 
tween them is much more important than the quantum 
corrections to the semiclassical limit. 

Effective dipole moment of the impurity and destruc- 
tion of the Neel order at 2% doping. It is convenient 
to rewrite Eq. © using the notation of the non-linear o~- 
model. Far from the impurity core, r> 1/k, the solution 



(e ■ r ) 

Sn = m6 = MM - v ' 



2nr 2 



M = 



y/2Zt 



8.2 
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(11) 



where n = 8n + no is the unit vector of antiferromag- 
netism, no = n(r = oo), m = [no X m], and m is the di- 
rector of the impurity. Here M is the effective dipole mo- 
ment of the impurity. Note that it is very large, M»l. 

The idea of destruction of the Neel order by ran- 
domly quenched dipoles was put forward by Glazman 
and Ioselevich |2,"J |. Detailed renormalization group cal- 
culations based on this picture have been performed by 
Cherepanov et al |26| and we use their results. In partic- 
ular an analysis of the experimental data by Keimer et 
al [2j for the in-plane correlation length at x < 0.02 was 
performed in |26|. This analysis shows that in order to 
explain the data and hence the destruction of the Neel 
order at x « 0.02 one needs to have a value of M which 
satisfies the following condition 



A = ^ = 20(1 ±0.3) 



(12) 



Here d = 2 is the dimensionality of the problem and M is 
the dimensionality of the vector m. In our theory Af = 2 
because m is orthogonal to no. Hence we conclude from 
JT2J that M exp = 8.9(1 ±0.15). This agrees well with the 
theoretical value (fill) . 

Structure of the insulating (spin glass) region and tran- 
sition into the metallic phase. Here we consider the range 
of doping 0.02 < x < 0.055 where the insulating spin 
glass state is realized. Since elastic incommensurate neu- 
tron peaks have been observed in this regime 0, 0, ' 
there are two characteristic length scales: cx 1/x, re- 
lated to the incommensurability, and the magnetic cor- 
relation length Im > h, related to the spin glass disorder 
(randomness) and reflected in the (inverse) width of the 
elastic neutron peaks. 

It is clear that in order to minimize the dipole-dipolc 
interaction energy at finite impurity concentration (and 
at zero temperature), the dipoles l|ll|l will align in such 
a way that all vectors e and m are the same. Such an 
alignment is possible in spite of the random positions and 
generates an average spiral 01 . Certainly around each 
dipole there are deviations from the average described 
by One can consider the average spiral as a self- 

consistent field created by all dipoles. To find the average 
pitch of the spiral let us consider a single dipole with field 
Sn given by (|TT|> in a background field n;, = n + Sni, 
where 

Sn b = -Xm b (e b ■ r) . (13) 

Here Sni, is the self-consistent field of the dipoles, e& = 
(1/V2, ±l/\/2) is a unit vector orthogonal to the face 
of MBZ, and A is a parameter. The interaction of the 
dipole with the background field is given by: 4f J (V<5n + 



\7Sn h ) 2 d 2 r 



^f(V5n) 2 d 2 ° 



J (VSnb) 2 d 2 r, which 
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simply amounts to p s J(V(5n)(V<5n;,)<i 2 r = — Mp s (m ■ 
mft)(e ■ e^). Clearly the interaction energy has a mini- 
mum at m = rrTf, and e = e;,. The total energy at a finite 
concentration x consists of the energy of each particu- 
lar impurity l|l(Jf) . the interaction energy, and the elastic 
energy of the background: 

E x = Ex~ p s CM\x + ^A 2 . (14) 

In the interaction energy term we have introduced the 
finite-size correction constant C. Indeed, Eq. Hll|) is 
valid only at very large distances from the impurity. 
However at a finite distance the effective dipole mo- 
ment is reduced, according to Eq. @, by the amount 
C = [l - e- 2Kr (l + 2/cr)] . Substituting r = 1/s/ttx and 
k w 0.4, we find for x = 0.03 - 0.05 the value C w 0.7. 
Minimizing l|14|l with respect to A we find A = CMx. 
Hence the average pitch is 

Zt 

SQ = \e b = C— x (1,±1) . (15) 

This expression determines the incommensurate shift of 
the neutron peak and agrees well with the experimental 
data flj since CZt/p s w 4.1. 

The last question we want to discuss is the micro- 
scopic origin of the correlation length Im- First we notice 
that without randomness we would automatically have 
Im =oo, whereas experimentally this quantity is about 
l M « 25-40A 0. It has been suggested in Ref. [p| that 
topological defects related to the random positions of im- 
purities can destroy the long-range spiral order. This is 
a possible scenario, however we suggest a different mech- 
anism. In our opinion even randomly distributed dipoles 
(fTTfl would create a true long-range spiral order, simi- 
larly to a 2D ferroelectric. However, our main observa- 
tion is that the situation is not fully described by the 
point-like dipoles. Each impurity has a finite size core 
(see Eq. ©) with diameter 1/k ~ 3 — 5 lattice spac- 
ings (depending on doping). Therefore, given a random 
distribution of positions, there is always a finite prob- 
ability of impurity overlap. As soon as the impurities 
overlap, a two-hole "molecule" is formed and the situa- 
tion changes dramatically. In the "molecule" the Pauli 
blocking starts to play a role and in order to minimize 
the energy the holes prefer to occupy different pockets in 
momentum space. If two pockets are occupied the (1,0) 
spiral has lower energy, and this is exactly what happen 
in the conducting phase 0. Hence such a "molecule" 
has a local spiral along (1,0) or (0,1) direction. This spi- 
ral frustrates the (1,1) background and there is always 
a finite concentration of such frustrating dipoles. Hence 
the "molecule" dipoles destroy the (1,1) background sim- 
ilarly to the way the "atomic" (single) dipoles destroy the 
Neel background at x < 0.02. One can consider these 
"molecules" as a precursor to the transition to the con- 
ductor where the (1,0) spiral is realized. According to 



this picture the spin-glass correlation length Im is large 
(but always finite) at very small x and it should decrease 
dramatically towards the percolation point x w 0.055 
where the "molecular" configurations are becoming more 
important. This is exactly what is observed in experi- 
ment, as seen in Fig. 6 of Ref. 0. In the superconduct- 
ing phase the magnetic correlation length should increase 
very rapidly, since theoretically it is infinity in the fully 
uniform, metallic phase Indeed, experimentally the 
correlation length quickly approaches the uniform limit 
7] (it is > 200A, x = 0.12). 

In conclusion, we have developed a description of 
the magnetic properties of undcrdoped I^-xSr^CuO.j, 
based on the extended t— J model. The theory describes 
the incommensurate elastic neutron scattering above and 
below the metal-insulator transition at x — 0.055. In par- 
ticular it explains why the incommensurate peak position 
rotates by 45° exactly at the insulator-metal transition. 
The theory does not contain any fitting parameters, and 
the positions of the neutron peaks both in the conducting 
J5J and in the insulating 1)15(1 phases, as well as the crit- 
ical concentration for destruction of the Neel order, fol- 
low from the calculated parameters of the extended t — J 
model. We also note that in L^-zSr^CuC^ static charge 
modulation (stripes) has not been directly observed, sug- 
gesting that it is very weak or not present at all. We thus 
believe that a theory based on spiral magnetic correla- 
tions and no charge order is fully sufficient to describe 
the phenomena in this material. 
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